home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / histogram / oper2d.c < prev    next >
Encoding:
C/C++ Source or Header  |  2001-11-01  |  4.3 KB  |  204 lines
  1. /* gsl_histogram2d_oper.c
  2.  * Copyright (C) 2000  Simone Piccardi
  3.  *
  4.  * This library is free software; you can redistribute it and/or
  5.  * modify it under the terms of the GNU General Public License as
  6.  * published by the Free Software Foundation; either version 2 of the
  7.  * License, or (at your option) any later version.
  8.  *
  9.  * This program is distributed in the hope that it will be useful,
  10.  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  11.  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  12.  * General Public License for more details.
  13.  *
  14.  * You should have received a copy of the GNU General Public
  15.  * License along with this library; if not, write to the
  16.  * Free Software Foundation, Inc., 59 Temple Place - Suite 330,
  17.  * Boston, MA 02111-1307, USA.
  18.  */
  19. /***************************************************************
  20.  *
  21.  * File gsl_histogram2d_oper.c: 
  22.  * Routine to make operation on 2D histograms. 
  23.  * Need GSL library and header.
  24.  * Contains the routines:
  25.  * gsl_histogram2d_same_binning check if two histograms have the same binning 
  26.  * gsl_histogram2d_add          add two histogram
  27.  * gsl_histogram2d_sub          subctract two histogram
  28.  * gsl_histogram2d_mult         multiply two histogram
  29.  * gsl_histogram2d_div          divide two histogram
  30.  * gsl_histogram2d_scale        scale histogram contents
  31.  *
  32.  * Author: S. Piccardi
  33.  * Jan. 2000
  34.  *
  35.  ***************************************************************/
  36. #include <config.h>
  37. #include <stdlib.h>
  38. #include <gsl/gsl_errno.h>
  39. #include <gsl/gsl_histogram2d.h>
  40.  
  41. /* 
  42.  * gsl_histogram2d_same_binning:
  43.  * control if two histogram have the
  44.  * same binning
  45.  */
  46. int
  47. gsl_histogram2d_equal_bins_p (const gsl_histogram2d * h1,
  48.                   const gsl_histogram2d * h2)
  49. {
  50.  
  51.   if ((h1->nx != h2->nx) || (h1->ny != h2->ny))
  52.     {
  53.       return 0;
  54.     }
  55.   {
  56.     size_t i;
  57.     /* init ranges */
  58.     for (i = 0; i <= (h1->nx); i++)
  59.       {
  60.     if (h1->xrange[i] != h2->xrange[i])
  61.       {
  62.         return 0;
  63.       }
  64.       }
  65.     for (i = 0; i <= (h1->ny); i++)
  66.       {
  67.     if (h1->yrange[i] != h2->yrange[i])
  68.       {
  69.         return 0;
  70.       }
  71.       }
  72.   }
  73.   return 1;
  74. }
  75.  
  76. /* 
  77.  * gsl_histogram2d_add:
  78.  * add two histogram
  79.  */
  80.  
  81. int 
  82. gsl_histogram2d_add (gsl_histogram2d * h1, const gsl_histogram2d * h2)
  83. {
  84.   size_t i;
  85.  
  86.   if (!gsl_histogram2d_equal_bins_p (h1, h2))
  87.     {
  88.       GSL_ERROR ("histograms have different binning", GSL_EINVAL);
  89.     }
  90.  
  91.   for (i = 0; i < (h1->nx) * (h1->ny); i++)
  92.     {
  93.       h1->bin[i] += h2->bin[i];
  94.     }
  95.  
  96.   return GSL_SUCCESS;
  97. }
  98.  
  99. /* 
  100.  * gsl_histogram2d_sub:
  101.  * subtract two histogram
  102.  */
  103.  
  104. int 
  105. gsl_histogram2d_sub (gsl_histogram2d * h1, const gsl_histogram2d * h2)
  106. {
  107.   size_t i;
  108.  
  109.   if (!gsl_histogram2d_equal_bins_p (h1, h2))
  110.     {
  111.       GSL_ERROR ("histograms have different binning", GSL_EINVAL);
  112.     }
  113.  
  114.   for (i = 0; i < (h1->nx) * (h1->ny); i++)
  115.     {
  116.       h1->bin[i] -= h2->bin[i];
  117.     }
  118.  
  119.   return GSL_SUCCESS;
  120. }
  121.  
  122. /* 
  123.  * gsl_histogram2d_mult:
  124.  * multiply two histogram
  125.  */
  126.  
  127. int 
  128. gsl_histogram2d_mul (gsl_histogram2d * h1, const gsl_histogram2d * h2)
  129. {
  130.   size_t i;
  131.  
  132.   if (!gsl_histogram2d_equal_bins_p (h1, h2))
  133.     {
  134.       GSL_ERROR ("histograms have different binning", GSL_EINVAL);
  135.     }
  136.  
  137.   for (i = 0; i < (h1->nx) * (h1->ny); i++)
  138.     {
  139.       h1->bin[i] *= h2->bin[i];
  140.     }
  141.  
  142.   return GSL_SUCCESS;
  143. }
  144.  
  145. /* 
  146.  * gsl_histogram2d_div:
  147.  * divide two histogram
  148.  */
  149.  
  150. int 
  151. gsl_histogram2d_div (gsl_histogram2d * h1, const gsl_histogram2d * h2)
  152. {
  153.   size_t i;
  154.  
  155.   if (!gsl_histogram2d_equal_bins_p (h1, h2))
  156.     {
  157.       GSL_ERROR ("histograms have different binning", GSL_EINVAL);
  158.     }
  159.  
  160.   for (i = 0; i < (h1->nx) * (h1->ny); i++)
  161.     {
  162.       h1->bin[i] /= h2->bin[i];
  163.     }
  164.  
  165.   return GSL_SUCCESS;
  166. }
  167.  
  168. /* 
  169.  * gsl_histogram2d_scale:
  170.  * scale a histogram by a numeric factor
  171.  */
  172.  
  173. int 
  174. gsl_histogram2d_scale (gsl_histogram2d * h, double scale)
  175. {
  176.   size_t i;
  177.  
  178.   for (i = 0; i < (h->nx) * (h->ny); i++)
  179.     {
  180.       h->bin[i] *= scale;
  181.     }
  182.   
  183.   return GSL_SUCCESS;
  184. }
  185.  
  186. /* 
  187.  * gsl_histogram2d_shift:
  188.  * shift a histogram by a numeric offset
  189.  */
  190.  
  191. int 
  192. gsl_histogram2d_shift (gsl_histogram2d * h, double shift)
  193. {
  194.   size_t i;
  195.  
  196.   for (i = 0; i < (h->nx) * (h->ny); i++)
  197.     {
  198.       h->bin[i] += shift;
  199.     }
  200.   
  201.   return GSL_SUCCESS;
  202. }
  203.  
  204.